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Abstract 

We present a theoretical discussion of the reversible parking prob- 
lem, which appears to be one of the simplest systems exhibiting glassy 
behavior. The existence of slow relaxation, nontrivial fluctuations, 
and an annealing effect can all be understood by recognizing that two 
different time scales are present in the problem. One of these scales 
corresponds to the fast filling of existing voids, the other is associated 
with collective processes that overcome partial ergodicity breaking. 
The results of the theory are in a good agreement with simulation 



data; they provide a simple qualitative picture for understanding re- 
cent granular compaction experiments and other glassy systems. 

PACS numbers: 68.45. Da, 61.43.-j, 64.70.Pf 
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1 Introduction 



Constrained dynamics of complex systems has long been a subject of exten- 
sive experimental and theoretical research. Certain important features, such 
as slowing kinetics, non-exponential relaxations, and memory effects are be- 
lieved to be generic for a wide class of systems ranging from structural and 
spin glasses to granular materials and traffic flows 0, [[3]. In spite of 
significant progress in this field there is still no general framework for the 
description of jamming and glassy phenomena. The development of a clear 
qualitative picture of them has been frustrated by the relatively high com- 
plexity of the considered systems. One could hope to boost the conceptual 
progress in this field by analyzing simple models capable of capturing the 
important features of glasses. 

In this paper we present a theoretical discussion of one of the simplest 
systems that exhibits glassy-like relaxation dynamics and a non-trivial fluc- 
tuation spectrum. This system is known as the Parking Lot Model (PLM) 
or the continuous Random Adsorption Problem, and it is defined as 
follows. Identical, unit length particles (cars) can adsorb on a line (curb) 
at rate k + per unit curb length. They can also leave the line with rate 
The desorption process is unrestricted while the adsorption is subject to free 
volume constraints, i.e. two cars cannot overlap (see Figure 1). This model 
can be applied in a straightforward way to random physical adsorption of 
large molecules. In addition, the PLM appears to be one of the most suc- 
cessful models for the description of density relaxation and fluctuations in a 
vibrated granular material. The possible reason for this is that the dynamics 
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of the PLM drastically depends on the available free volume, just as in the 
case of granular materials or structural glasses. 

The dynamics of the original version of the PLM, in which the particles 
adsorb irreversibly (i.e. k_ = 0), has been well understood as long ago as 
1958 by Renyi [f|. He found that the system jams at density p c ~ 0.75, and 
the way it approaches this state is given by the following formula: 



The late-stage asymptotics of this result is the power law relaxation: 



The desorption process introduced by Krapivsky and Ben-Nairn ||, re- 
sults in even richer physics. In a recent paper on granular compaction, Nowak 
et. al. presented the results of simulations on the PLM. These simula- 
tions and the experimental granular compaction data have many important 
features in common. In particular, the average coverage of the curb as a func- 
tion of time, shown in Figure 2, is very similar to the density relaxation curve 
for the vibrated sand. Once in the steady state, the finite size of the system 
results in considerable density fluctuations. Insight into their dynamics is 
provided by the power spectrum as seen in Figure 3. This figure illuminates 
one of the most remarkable properties of the parking lot model, i.e. that it 
exhibits two very different time scales at high k values (high density). These 
time scales appear in the power spectrum as two corner frequencies, one 
at high frequency and one at low. The low frequency corner is Lorentzian, 
which indicates that it can be associated with exponential relaxation at a 




(1) 



P(t) ~ Pc ~ 1 



It. 
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single time scale. The high frequency "corner", however, shows an unusual 
"hump" that indicates that the density relaxation cannot be described with 
only one time scale. It is these features in the power spectrum that led 
Nowak et. al. to present the parking lot model in conjunction with the sand 
experiment, which shows a similar non-trivial fluctuation spectrum. The ex- 
istence of several relaxation time scales is a signature of the partial ergodicity 
breaking exhibited by the model, i.e. its high-frequency evolution does not 
allow the system to explore all the configurational space. In this sense the 
observed behavior of the PLM may be relevant not only for understanding 
the particular granular compaction experiment but also for the whole class of 
systems exhibiting glassy relaxation dynamics. Below we focus on developing 
an analytical description of the PLM capable of capturing these intriguing 
features and revealing the underlying physics. 

2 The Limitations of the Mean Field Ap- 
proach 

In their original work on the reversible PLM ||, Krapivsky and Ben-Nairn 
proposed a mean-field description of the problem that can be essentially 
expressed in terms of the following master equation for the average density, 
P- 
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Here the first term represents the adsorption rate and takes into account its 
strong dependence on the available free volume, the second term corresponds 
to unrestricted desorption. The above equation indeed captures some of 
the features of the PLM. In particular, it results in an equation for the 
equilibrium density p eg , that is consistent with the numerical data: 



Although this mean-field approach does result in typical jamming dynam- 
ics, the predicted density relaxation rate is several orders of magnitude faster 
than that observed in the simulations 0. Moreover, the mean field descrip- 
tion fails to capture the most interesting feature of the PLM, its fluctuation 



different from a regular linear response equation for the density near equilib- 
rium, i.e. it should have a single relaxation time for small fluctuations. 

An essential feature that the above mean-field approach overlooks is the 
strong correlation between adsorption and desorption events. This correla- 
tion is familiar to anyone who has ever attempted to park in a big city. Cars 
do not leave the curb very often but as soon as they do an incoming car 
rushes in to take the newly created space. This paired adsorption-desorption 
process does not change the density of the system; it is equivalent to merely 
sliding a car in its parking space. One can explicitly emphasize this by the 
introduction of a slow (adiabatic) variable Z = N + iV*, which is the sum 
of the number of adsorbed cars, N, and the number of voids large enough 
to fit at least one particle, N*. If the system is sufficiently dense, most indi- 






cannot result in anything 
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vidual adsorption or desorption events do not change the parameter Z. The 
existence of this slow variable results in the observed separation of the relax- 
ation time scales. The fast modes of the density relaxation correspond to the 
evolution with nearly constant Z, while the slow dynamics is determined by 
low-probability events, which result in a drift of Z . In other words, on short 
time scales the ergodicity of the PLM is broken: the system can explore only 
the part of the configurational space corresponding to a constant Z . 

The fast and slow modes are expressed as two typical time scales in the 
density fluctuations and also lead to a two-stage relaxation process of the 
density evolution from p = to p eq . In the case of a very small desorption 
rate k- the system follows the universal irreversible-adsorption curve, Eq. 
([]]) until it jams at p c ~ 0.75. Afterward, it slowly evolves towards the 
equilibrium density. By using the above result for 

Peq cL function of k = 

k + /k^ and knowing p c one can construct the "kinetic phase diagram" of the 
system (see Figure 4): if p eq (k) < p c the system reaches the equilibrium 
density before it jams, otherwise it relaxes to p c and then slowly creeps 
toward equilibrium. It would be tempting to associate the critical point 
k c ~ 60, which separates the two regimes, with a glass transition. However, 
this term is traditionally reserved for the hypothetical point at which the 
kinetic coefficients of a system would go to zero (in our simple case this is 
the point k = oo). Thus, the zone of partial ergodicity breaking appearing 
is an analog of what is conventionally called a supercooled liquid 
rather than a glass. 

An important implication of the above picture of the PLM kinetics is the 
existence of an annealing effect, which is typical for classical glassy systems 
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and has also been observed in the granular compaction experiments || . The 
idea is that one can overcome the slow kinetics of the supercooled system by 
"heating" it (by decreasing k below k c ), then slowly cooling it (increasing k), 
so that the system would follow the reversible equilibrium curve p eq (k). 



3 The Fast Dynamics of the System 

As we have pointed out above, the fast dynamics of the PLM is dominated 
by a two part process: a desorption that takes place at rate leaving a void 
and an adsorption that occurs at rate k + z, where z is the void size minus 
a car length. This two-part process is equivalent to simple replacement of 
one adsorbed particle with another or a sliding of a single particle in its 
own space. As a result, the separations between a given particle and its two 
nearest neighbors, x\ and 22, change randomly to the new values x[ and x' 2 
in such a way that their sum remains constant: x\ + X2 = x\ + x' 2 = z. This 
sliding process creates the necessary mixing to ensure that the steady-state 
distribution of the spaces between cars, p(x), corresponds to the maximal 
entropy for a fixed sum of all the separations. Thus, the distribution function 
can be obtained by maximizing the functional 

&[p] = ~~ J dx (p(x)logp(x) + p(x)^\ (4) 

with respect to p(x), where 1/5 is the Lagrange multiplier conjugate to x. 
Naturally, the resulting formula for p(x) is similar to the Boltzmann distri- 
bution: 

p(x) = I exp ( (5) 



6 r V 5 



The parameter 5 is equal to the average spacing between cars, i.e. 

, l-p 



(6) 



P 

We now proceed with the calculation of the high frequency part of the fluc- 
tuation spectrum. The coupled adsorption-desorption process can be viewed 
as a relaxation-excitation of the conventional two-state (telegraph) system, 
which can switch from one state, '0', to another one, T', with a characteristic 
time t and then relax back to '0' with the characteristic time r. The power 
spectrum of telegraph noise has been investigated by Machlup |j7J and is 
given by 

S(u) = — — ^ r—^ 1 — x (7) 

where v — 1/t + 1/t. 

The fast dynamics of the density is the superposition of these nearly 
independent telegraph modes, whose number is essentially the number of cars 
(or, more precisely, Z). An important feature of the adsorption-desorption 
modes is that while their "excitation rate", 1/t = is uniform over the 
system, the relaxation time is the property of an individual excited state, i.e. 
the adsorption rate is proportional to the free volume available at a given 
void: 1/r = k + (xi + x 2 ). As a result, the fluctuation spectrum deviates from 
the simple Lorentzian form, Eq (^): 

Sh(u) = {pup-u} = - 2 = — / 8 

Here L is the total length of the system. Note that here and below we assume 
that k = k+/k-. 3> 1, so that the typical adsorption process is much faster 
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than the desorption. As a result, the desorption rate determines the ampli- 
tude in the above expression, while the (non-uniform) adsorption dominates 
its frequency dependence (i.e. v ~ 1/r = k + (xi +x 2 )). Thus, the calculation 
of the power spectrum for the parking lot model reduces to finding the distri- 
bution of relaxation times for adsorption, and then weighting the spectrum 
for telegraph noise with this distribution. As shown above, the distribution 
function for the inter-car spaces is exponential. Assuming that the separa- 
tions between a given car and its two neighbors, x\ and X2, are uncorrelated, 
we obtain the following distribution of relaxation rates, v = k+(x\ + x 2 ): 



is the characteristic frequency of the fast relaxation. 

One can now substitute the distribution function for v into the above 
expression for the high-frequency fluctuation spectrum, Eq. (|8]). The result 
of the numerical integration of this formula over v is in excellent agreement 
with the simulations, see Figure 3. We now have an explanation for the hump 
in the spectrum at high frequency; it is caused by a relatively broad distri- 
bution of relaxation times for the adsorption process. Further examination 
of our result illuminates the connection between the PLM and other, more 
complicated systems exhibiting glassy dynamics. The above distribution of 
relaxation rates does not just distort the Lorentzian form near the character- 




(9) 



Here 



oo H = k + 5 



(10) 



10 



istic frequency, ujh = k + 5, but rather affects the spectrum at all frequencies 
below ujh- The asymptotic behavior of the spectrum at u <C ujh is 

S(u)~\ogu (11) 

Hence, the spectrum never recovers the Lorentzian-like plateau regime at 
low frequencies. The logarithmic behavior is reminiscent of the power-law 
spectrum typical for glasses at mesoscopic frequencies (it corresponds to f3- 
relaxation, [3J). Because the logarithmic behavior is a result of the prob- 
lem being one-dimensional, one could expect real power-law behavior in the 
spectrum in higher dimensions. A simple assumption that the adsorption 
probability is proportional to the free volume associated with the newly cre- 
ated void, combined with the natural exponential distribution for the free 
volume, would result in the 1/uj spectrum at mesoscopic frequencies: 

r°° exp(-u/u )du 
S h (uj) ~ / =- — - 1/u , 10 < u (12) 

JO UO + V 

Note that in the one-dimensional case discussed above, the free volume 
is a sum of two presumably independent variables, x\ and 22, each of which 
has an exponential distribution. As a result, the distribution function for the 
free volume vanishes near zero, i.e. the probability of finding a long-living 
excited state is strongly suppressed in the one-dimensional case compared to 
the higher dimensions. 
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4 The Slow Dynamics of the System 



In the previous section we have discussed the fluctuations of the density on 
short time scales, over which the system remains ergodically broken. Now we 
proceed to a discussion of the low-frequency part of the fluctuation spectrum 
associated with the change of the slow parameter Z. This change is caused 
by collective events; the rearrangement of a state corresponding to a given Z 
is dominated by a two-car process. In one process, responsible for decreasing 
the "ground state" Z by 1, two adjacent cars leave and a single one (a 
"bad" parker) comes in their stead, hogging the space (see Figure 5a); the 
opposite process results in adding an extra car to the lot: a car exits and 
leaves a large space-big enough for two cars, provided that the new cars are 
"good" parkers (figure 5b). How do these collective modes affect the power 
spectrum? We now calculate the rates of the above two-car processes. The 
rate of the " — 1" process has three contributions. First, a car must leave, and 
the corresponding "trial rate" is just A;_ per particle. Then, an adjacent car 
must leave before the hole left by the first car fills, which gives a "waiting" 
factor 2/c_ / °° J °° ex"p{—vt)dtj{y)dv. Finally, the big hole must be blocked 
by a "bad" parker (up to a correction of order of 5, the probability of this is 
unity). Thus the overall rate is 

2k 2 , . 

= M (13) 

The opposite, process, has the same trial frequency, k_ per car. 

The void left by the car (its length is xi + x 2 + 1) must be large enough for 
two cars, which gives the factor 
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1 + A 

f o eX P 



1 + A 



f . jrfA (14) 

Note that the first incoming car must park with precision A = x\ + X2 — 1 
in order to leave enough space for the second car. The probability of this 
happening contributes a factor of 2A/(1 + A) to the overall "+l"-rate: 





r i + ai 


la ^ 6XP 
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dA = 2k_exp(-l/5) (15) 



The density changes with time according to the following Equation: 

p = - v- x )p + Tj(t) (16) 

Here rj is the noise originating from the fact that the density changes by 
discrete one-particle steps. (77 (t) ) = 0, because the average evolution is 
given by the interplay of the and "—1" kinetic terms in the above 

equation. Since adding an extra car or removing one at the moment of time 
to corresponds to N = ±5(t — to), and since there is no obvious mechanism 
for the correlations between such processes, 

Using Eq. ( Jl6| ) , we find that the equilibrium density is determined by the 
condition: 

= 5 eq ex P (-l/5 eq ) (18) 

where 5 eq = p~ q l — 1. This result coincides with the mean field one, Eq. 
(|3p. We note that the mean field approach ignores the adsorption-desorption 
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correlations; this would be a reasonable assumption for the model with strong 
diffusion of the adsorbed cars. Since the diffusion cannot shift the equilibrium 
properties, it is not surprising that the mean field approach gives the correct 
value of p eq . As to the description of the PLM kinetics, the rates v_\ and z/ +1 
in Eq (|T6|) differ by an exponentially small factor, 2exp(— l/5 eq ) from their 
mean-field analogs. The above equation describes only the slow evolution of 
the jammed state toward equilibrium. Thus, the overall density relaxation 
curve of the "cold" system (at k > k c ) consists of the classical fast regime, Eq. 
(HD, resulting in a jamming at p = p c ~ 0.75, and the desorption-promoted 
final stage, discussed here. Such a combination of the two theoretical results 
is in agreement with the simulation data, as is shown in Figure 2. The same 
figure shows that the mean field curve does not capture the two stage nature 
of the relaxation dynamics and is inadequate for the description of the cold 
system (k = 10 4 ). However, the mean field may be used for the description 
of the single-stage relaxation of the "hot" system. 

By expanding Eq. (|16|) near the equilibrium density, one can determine 
the relaxation frequency of the system, uj^ and the spectrum of the low- 
frequency fluctuations: 

2k 

u L = -±exp(-2/6 eq ) (19) 

Oeq 

= 2k -A,«P(-i/f ,) (20) 

7iL (or + ujjj 

By combining this expression with the earlier result for the high-frequency 
fluctuations, one obtains the following analytic form for the entire power 
spectrum of the PLM: 
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S(u) = S H (v) + S L (cu) 



irL 



-Peg 



{I 



■oc 



exp(— v/ujh) vdv 2exp(— 1/S eq ) 
UJ 2 + v 2 oj 2 h UJ 2 + uj\ 




This result agrees amazingly well with the simulation data as shown in 
Figure 3. Figure 6 shows how the characteristic frequencies ojh and uj l 
depend on the control parameter k = k + /k_ (or, equivalently on p eq ); the 
theoretical calculation is again in good agreement with simulation. 

5 Conclusions 

We have presented a theoretical discussion of the Parking Lot Model, which 
appears to be a very simple glassy system, perhaps the simplest. We have 
identified two time scales in the problem: one associated with a simple re- 
laxation of voids and the other corresponding to the collective (two-particle) 
processes responsible for the rearrangement of the "ground state" (the state 
that the system can reach by an instant filling of all currently available 
voids). In the limit of weak desorption, corresponding to a large difference 
between the two times, the relaxation of the system toward its equilibrium 
density is a two-stage process: first, it reaches the universal jamming density 
p = p c ~ 0.75 as if there were no desorption at all, then it slowly relaxes 
to Peq(k) via collective rearrangements. This two-stage relaxation feature 
disappears in the regime of strong desorption (k < k c ), and we identify the 
crossover point k c with p eq (k) = p c . 

At k > k c , on times shorter than the longest characteristic scale {ojj}) the 
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evolution of the system is non-ergodic; this regime is analogous to a super- 
cooled liquid. The system evolves by jumping between the metastable states 
corresponding to different values of the parameter Z. Note that the lifetime 
of these states does not grow with the system size, but rather decreases. 
Since the probability of the rearrangement of the ground state is v + \ + v_\ 
per particle per unit time, its lifetime is inversely proportional to the number 
of cars iV ~ Z: r^ 1 = Z(u +1 + V-\). In this sense, the free energy landscape 
of PLM is similar to that of structural glasses. Another similarity between 
the PLM and glassy systems is the possibility of accelerating relaxation by 
means of annealing. 

The existence of two characteristic time scales is responsible for the in- 
triguing form of the fluctuation spectrum of the reversible parking problem. 
The slow fluctuations are described by a single Lorentzian with the relax- 
ation frequency ujl associated with two-particle rearrangements. The fast 
dynamics is a superposition of many single-particle adsorption-desorption 
modes. In a sense, the non-Lorentzian form of the high-frquency part of the 
power spectrum is a reflection of the deviations of the local density from its 
average value, i.e. it is a signature of long-living disorder. In this form, our 
observation may be relevant for understanding the non-trivial behavior of the 
relaxation spectrum of other glassy systems on mesoscopic frequencies. The 
distribution of the relaxation rates of the single-particle excitations results 
not only in the distortion of the Lorentzian in the vicinity of the character- 
istic corner frequency uh, but also in an interesting logarithmic behavior of 
the power spectrum at lower frequencies (S ~ logo;). We have suggested 
that this feature is reminiscent of the power law spectrum corresponding to 
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/^-relaxation in classical glasses and that such a power-law behavior could 
be reproducible in PLM at higher dimensions. 
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Figure Captions 

Figure 1. The parking lot model. 

Figure 2. Simulation data (circles) and theoretical results for density evo- 
lution in the "cold" (k = 10 4 ) amd "hot" (A; = 10) regimes. Note that 
the mean field (dashed line) is quite adequate for the description of 
the "hot" system but it fails to describe the slow dynamics at the "su- 
percooled" regime for k = 10 4 . The solid line is the combination of 
the irreversible parking curve, Eq (1), describing the fast stage of the 
relaxation (for k = 10 4 ), and our result, beginning at t = 100, for the 
later slow dynamics. 

Figure 3. The power spectrum of density fluctuations near the equilibrium 
for k = 10 3 and k = 10 4 . Circles and solid lines represent the simulation 
data and our analytic results, respectively. 

Figure 4. "Kinetic Phase Diagram" of the system. The solid line corre- 
sponds to the reversible equilibrium line p eq (k), the dashed line shows 
the jamming density above which the ergodicity is partially broken. 

Figure 5. The slow collective rearrangements corresponding to (a) adding 
and (b) subtracting one car. 

Figure 6. The low and high characteristic frequencies, u>l and u H , as func- 
tions of the equilibrium density. The dashed and solid lines represent 
the theoretical results, the squares and the circles are the simulation 
data for ul and uh, respectively. 
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